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We study renormalized quenched strong-coupling QED in four dimensions in ar- 
bitrary covariant gauge, in the Dyson-Schwinger equation formalism. Above the 
chiral critical coupling, we show that there is no finite chiral limit. This behaviour 
is found to be independent of the detailed choice of proper vertex, provided that 
the vertex is consistent with the Ward-Takahashi identity and multiplicative renor- 
malizability. The finite solutions previously reported lie in an unphysical regime 
of the theory with multiple solutions and oscillating mass functions. This study is 
consistent with the assertion that strong coupling QED4 does not have a continuum 
limit in the conventional sense. 



1 Introduction 

The abelian nature of quantum electrodynamics (QED) in many ways makes 
it a much simpler system to study than a nonabelian theory such as quan- 
tum chromodynamics (QCD). For this reason it has been the subject of many 
nonperturbative studies, which hav^aa-their long-term goal a detailed under- 



standing of nonperturbative QCD. Elutil On the other hand, strong-coupling 
QED4 is widely anticipated to behave unconventionally in the continuum limit 
and for this reason is^theory of considerable interest in its own right. 

In previous workHEl we introduced a numerical renormalization procedure 
for the Dyson-Schwinger equations, and applied it to QED4 with a quenched 
photon propagator, under momentum cutoff regularization. The central result 
of those works was to demonstrate that the numerical renormalization proce- 
dure works extremely well and allows the continuum limit (cutoff A — > 00) to 
be taken numerically, while giving rise to stable finite solutions for the renor- 
malized fermion propagator. 

Here we describe ..studies of the chiral limit in renormalized quenched 
strong-coupling QED 4 pl using a photon- fermion vertex that satisfies the Ward- 
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Takahashi Identity (WTI) and makes the fermion DSE multiplicatively renor- 
malizable. We find that for couplings above the chiral critical coupling, keeping 
the bare mass mo (A) = as the cutoff is relaxed results in either a dynamical 
mass function which is identically zero, or a dynamical mass function which di- 
verges with the cutoff. The previously reported finite solutions showed damped 
oscillations in the dynamical mass functions at large p 2 , which suggested that 
they were unphysical. Further, we show here that for a given supercritical cou- 
pling and the same bare mass mo (A), it is possible to have multiple solutions 
corresponding to different renormalized masses m(p). We conclude that above 
the chiral phase transition, quenched strong-coupling QED in four dimensions 
does not have a chiral limit in the conventional sense. 



2 Formalism 



We write the fermion propagator as 

o ( , _ Z(p 2 ) _ 1 
[P> i>~M{f) A(p 2 ) i>-B(p 2 ) K ' 

where we refer to A(p 2 ) = 1/Z(p 2 ) as the finite momentum-dependent fermion 
renormalization and where M(p 2 ) = B{p 2 ) / A(p 2 ) is the fermion mass func- 
tion. Dynamical chiral symmetry breaking (DCSB) occurs when the fermion 
propagator develops a nonzero scalar self-energy in the absence of an explicit 
chiral symmetry breaking (ECSB) fermion mass. We will refer to coupling 
constants strong enough to induce DCSB as supercritical and those weaker as 
subcritical. 

The DSE for the renormalized fermion propagator, in an arbitrary covari- 
ant gauge, is 

S~\p) = Z 2 Qi, A)[ff - m (A)] - iZ 1 ( f ,, A)e 2 J ^^S(k)T^k,p)D tlu (q) ; 

(2) 

here q = k — p is the photon momentum, ji is the renormalization point, and A 
is a regularizing parameter (taken here to be an ultraviolet momentum cutoff), 
and mo (A) is the regularization-parameter dependent bare mass. We work in 
the quenched approximation, i.e., the photon propagator is bare, so that for 
the coupling strength and gauge parameter we have a = e 2 /4ir = ay = e^/An 
and £ = Co- n| -. 

Among other general requirements EH3, the vertex Ansatz is constrained 
by the Ward- Takahashi Identity (WTI), which is necessary for gauge invari- 
ance and guarantees the equality of the propagator and vertex renormalization 
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constants, Z2 = Z\. Ball and ChiuH have described the most general fermion- 
photon vertex that satisfies the WTI; it consists of a longitudinally-constrained 
(i.e., "Ball-Chiu") part Fg C , which is a minimal solution of the WTI with no 
artificial kinematic singularities, and a set of eight vectors T^(k,p), which 
span the transverse subspace and are also kinemalically regular as k — > p (im- 
proved versions of the T/* can be found in Ref . □) . A general vertex is then 
written as Y^{k,p) = T^ c (k,p) + £),- =1 T i(k 2 ,P 2 , q 2 )Tl l (k,p) , where the n are 
functions that must be chosen to ensure the correct discrete symmetry prop- 
erties, perturbative limits etc For our numerical examples we have used the 
Curtis-Pennington vertexEJ, however our results and conclusions regarding the 
chiral limit will not depend on the specific vertex Ansatz as long as it respects 
multiplicative renormalizability. 

The renormalization procedure is entirely analogous to one-loop pertur- 
bative renormalization: One first determines a finite, regularized self-energy, 
£'(//, A; f>), depending on both the regularization parameter A and the renor- 
malization point //. The DSE for the renormalized fermion propagator can be 
written in terms of E'(/x, A;p) as 

S-'ip) = 220i,A)^-mo(A)]-E>,A;p) 

= ji- m (n)-X(n;p) = A(p 2 )0-B(p 2 ) 1 (3) 

where £(/it;p) denotes the renormalized self-energy. The self-energies are de- 
composed into Dirac and scalar parts, E'(/z,A;p) = T,' d (fi,A;p 2 ) ^+E' s (/x, A; p 2 ), 
(and similarly for E(/i,p)). By imposing the renormalization boundary condi- 
tion, 

S-^)\ p ^ =jf-m(A*), (4) 

one finds relations for the scalar and Dirac self-energies, propagator renormal- 
ization constant Ziiji, A), and bare mass. A mass renormalization constant 
can also be defined as Z m (fj,,A) = mo(A)/m(/x), i.e., as the ratio of the bare 
to renormalized mass. 

We approach the chiral limit by setting the bare mass to zero and removing 
the regularization, i.e., maintaining m (A) = while taking the lirnitJV — > 00; 
this is the same definition used in nonperturbative studies of QCEBEffi. Obvi- 
ously, any limiting procedure where we take too (A) — > sufficiently rapidly as 
A — » 00 will also lead to the chiral limit. 

The behavior of the solutions under a renormalization point transformation 
is governed by a "renormalization group flow" : under any transformation fi — > 
//, we must have for all p 2 

m(h';p 2 ) = m(h;p 2 ) = m(p 2 ) , 
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so that, for the fermion propagator, S (/x' ; p) / S (fx; p) — Zi(n, &) / Zi(n' , &) in 
the usual way. We have verified Eq. (|^) to within 1 part in 10 4 for our calcu- 
lations. 

The DSE is solved in Euclidean momentum-space, after separation into 
Dirac-odd and -even parts describing the finite renormalization A(p 2 ) and the 
scalar self-energy B(p ). We use a "gauge-covariance- improved" momentum 
cutoff scheme, which is described in detail in Appendix A of Rcf. u Choosing 
the renormalized mass, m(fi), and then solving for the bare mass, mo (A), leads 
to well-behaved finite solutions for ALp 2 ) and M(p 2 ), which do not vary as we 
take the continuum limit (A — > oo) I 



3 Results and Conclusions 



A(p ) for several UV cutoffs 
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M(p") for several UV cutoffs 
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Figure 1: The behaviour of the finite renormalization A(p 2 ) and the mass function M(p 2 ) as 
a function of the ultraviolet cut-off A for mo (A) = 0. These solutions were for renormalization 
point p? = 10 4 , coupling a = 1.15, and gauge parameter £ = 0.25. As A — > oo we find 
A(p 2 ) —> 1 for all p 2 and M(p 2 ) diverges with A. 



Fig. [j] shows the non-oscillating solution A(p 2 ) and M(p 2 ) for supercritical 
coupling and with mo (A) = 0, for a wide range of values of A. It is clear 
from these numerical solutions that in the continuum limit (i.e., A — > oo) and 
for supercritical coupling, we find A(p 2 ) — > 1 for all p 2 and a mass function 
M(p 2 ) which diverges proportionally to A. Conversely, for subcritical coupling 
(i.e., a < ct c ), the chiral limit does exist and we find simply that A(/i;p 2 ) = 
(p 2 /li 2 )~ a Z/ 47T . It seems clear from these numerical studies that above critical 
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coupling there is no finite chiral limit in the continuum quenched theory in any 
covariant gauge, even in the presence of the renormalization procedure. 

While these conclusions were based on a numerical study with a specific 
choice of vertex, one can easily construct a general argument which applies 
irrespective of this choice: Consider any vertex_T which satisfies the WTI and 
which leads to multiplicative renormalizabilityB. It automatically follows that 
M(p 2 ) and A(p 2 )/Z2{n, A) are renormalization point independent, as in Eq. 
(|J) . We can then define dimensionless quantities by appropriately scaling with 
A, i.e., n ee fj,/A, p 2 ee p 2 /A 2 , M(p 2 ) ee M(p 2 )/A, A{p 2 ) ^A(^p 2 )/Z 2 (^ A), 
and mo = mo(A)/A. Note that the renormalization condition A(/i;/i 2 ) = 1 
automatically determines i?2(/i,A) for a given solution for fixed A. The di- 
mensionless functions M(p 2 ) and A(p 2 ) can only depend on dimensionless pa- 
rameters, i.e., a, £, fi, and ffiQ. Furthermore, since for any fixed A they are 
independent of fj, (recall that we are working only in the quenched approxima- 
tion), then it follows that they must in turn be independent of p,. Now for any 
finite A and the choice mo (A) = we have m,o = 0. Hence solving for any /i 
and A with mo (A) = allows us to form M(p 2 ) and A(p 2 ), from which we can 
read off the solutions for A(p 2 ) and M(p 2 ) for any other \x and A with van- 
ishing bare mass using the above rules. We see then that the resulting M(p 2 ) 
must diverge with A. So if M(p 2 ) ^ 0, it must diverge with A as was found 
numerically. To summarize, above critical coupling any vertex which satisfies 
the WTI and leads to multiplicative renormalizability will lead to a diverging 
mass function in the chiral limit for quenched QED in four dimensions. 

It remains to discuss the meaning of the well-behaved finite solutions. Our 
conventional thinking about the chiral limit would imply that since ECSB 
should only increase the mass function above that found in its absence, then 
any solution above critical coupling which also has ECSB must also correspond 
to a divergent mass in the continuum limit. Thus the finite solutions do not cor- 
respond to the chiral limit nor to any conventional concept of ECSB. Further, 
the finite solutions show damped oscillations in the dynamical mass function, 
periodic in In (p 2 ), and a corresponding oscillatory behavior for the bare mass. 
Thus it is possible for a given set of parameters a, mo (A), fj,, £, and A to admit 
more than one finite solution. We explicitly show this in Figs. ^ and pL where 
for a given set a, mo (A), [i, £, and A there are distinct solutions. As A — > oo 
the number of simultaneous solutions becomes infinite. It seems likely that this 
behavior for the finite solutions is independent of the detailed vertex choice 
and furthermore that the same behavior can be induced in QCD by forcing 
the renormalized mass into an unphysical regime (by choosing m(/x) EE M([i 2 ) 
below the value corresponding to the chiral limit). It also seems likely that 
unquenching the theory will not remove this rather undesirable behaviour in 
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Bare vs Renormalized Mass 
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Figure 2: The relationship between the bare mass (mo (A)) and the renormalized mass 
(m(fM)) for the renormalized finite solutions. These result from solving for a given m(/i) 
and extracting mo (A). The parameters for these solutions were the renormalization point 
fi 2 = 10 4 , a = 1.25, gauge parameter £ = 0.25 and A 2 = 10 14 . The dashed horizontal line 
shows, e.g., that for mo(A) = 0.1 there are three solutions. (The dashed vertical lines are 
merely to guide the eye on this back-to-back log scale.) 

QED4, since the running coupling increases with scale rather than decreasing 
as in QCD. These latter conjectures are the subject of current investigation^. 
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